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ABSTRACT 

This paper develops a model for computing, the proba- 
bility of kill for an air defense complex composed of 
antiaircraft automatic weapons, radar controlled guns, and 
missile batteries. Two dimensional terrain was used to 
evaluate the model. The probabilities were determined at 
major terrain points along the route of approach to the 
vital area for altitudes of up to 3000 feet above terrain. 
The curves of probability of kill versus altitude were 
found to be dependent on terrain, air defense tactics, and 
weapon system parameters. A survivability index is calcu- 
lated by combining the probabilities of kill with a pilot 
visual navigational probability. The resulting curves of 
survivability index versus altitude were found to be non- 
linear requiring a nonlinear programing technique to solve 
for the altitude of optimal survivability index within 
aircraft flight path constraints. The nonlinear solution 
was not included in this work. 
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CHAPTER I 



During the first year of the United States' involvement 
with the war in South East Asia, air power was employed 
extensively against many targets in North Vietnam. The 
commanders of attack aviation units were required to plan 
and execute strikes against the same targets over and over 
again. Initially air defense of these targets was sparse 
to nonexistent. As time passed, air defense weapons were 
supplied from communist bloc countries until the problem 
presented to the attack force commander itfas what to do 
about higher and higher attrition rates per raid. Specific- 
ally, strike planners were seeking a solution to the problem 
of finding the route of approach and altitude to fly into 
the target which would minimize aircraft losses from air 
defense means. One method of finding an answer to the 
problem is to analyze the losses from previous raids. This 
method is costly in lives and equipment and may not have 
produced data describing results from many altitudes and 
routes into the target. Using previous experience alone, 
one may never discover the optimal route. A major factor i 
in the effectiveness of land air defenses is the terrain. 
Very often terrain provides a natural route of approach 
which in itself minimizes the capability of the air defense 
efforts. Sometimes this fact is obvious and in other 
situations it is not. It would be desirable in such a 
situation to provide the attack air force commanders with a 
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method that determines the optimal approach altitude and 
route that minimizes the air defense effort. The method 
would generate an actual flight profile for the pilot. As 
new data on the enemy’s air defenses are ascertained, up- 
dated flight profiles would be provided in the time period 
between strikes. The method referred to here is a digital 
computer model which evaluates all possible routes of 
approach into a target. For each route, the model evaluates 
the terrain, the air defenses situated therein, and the pilot 
navigational problem for all feasible altitudes up to a 
maximum value, and then determines the altitude along the 
route which maximizes aircraft survivability. With a 
functioning computer program which computes such a flight 
profile, new information could be inserted into the program 
and new results obtained in a matter of hours. Such a tool 
would be of tremendous value to attack air force commanders. 

The development of a mathematical model and computer 
solution for generating a flight profile which maximizes 
strike aircraft survivability was accomplished in a thesis by 
Lieutenant Colonel W. S. Miller, Jr., USMC, and Major E. E. 
Brown, USMC, reference 1. In this work a model for calcu- 
lating probabilities of detection of strike aircraft over 
two dimensional terrain was combined with calculated proba- 
bilities of detection of ground navigational targets to 
generate a functional relationship of cost versus altitude 
of the flight path above terrain. The term cost is defined 
to be a probability index which is a linear combination of 
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probability of radar detection and the probability that the 
ground navigational target was not detected at each alti- 
tude above terrain from 100 feet to 3 >000 feet. Using this 
function as an objective function and calculating aircraft 
flight path constraints for each terrain point, the problem 
of determining that altitude which minimized the cost was 
solved as a linear program. The probability of detection 
mentioned above was calculated by computing at each terrain 
point and altitude the ratio of terrain visable to the 
total terrain within the coordinate system. As altitude 
increased, less terrain was masked from the point in question. 
Therefore, the probability of detection was found to increase 
in approximately a linear fashion. This model assumed that 
a radar within the air defense system could be placed at 
any and all points along the terrain under the aircraft 
flight path. In other words, radar locations were considered 
to be uniformly distributed over the entire route of approach. 
To be sure, this method of calculating the probability of 
detection produces an indication of the effect of altitude 
on terrain masking in evaluating air defense capability. 

This paper extends the work of Miller and Brown. Since 
the uniform distribution of radar sites is not realistic and 
the method of calculating probability of radar detection is 
at best only an indication of air defense capability, a 
mathematical model and computer program will be developed to 
calculate the probability of kill for three types of air 
defense weapons sited in the terrain. Specifically, the 
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problem to be solved will be to emplace an air defense com- 
plex in a section of terrain in order to calculate at each 
terrain point and altitude above terrain from 100 feet to 
3,000 feet along the aircraft flight path the cumulative 
probability of kill for an aircraft flying straight and 
level up to each terrain point. The terrain points con- 
sidered are those within the maximum effective range of the 
longest ranged air defense unit to the defended area. The 
problem reduces to finding the length of course line exposed 
to radar detection and within the effective range of the 
weapon systems in order to determine the number of rounds 
that may be fired from which the probability of kill is 
determined. The model evaluates the terrain masking on 
detection range and firing time. The problem is developed 
for three dimensional terrain. The computer model was 
programmed for three dimensional terrain. No digitalized 
three dimensional terrain was found for use in this model. 
The computer program was accordingly verified using only 
two dimensional terrain which was available. That portion 
of the program which evaluates the three dimensional terrain 
was not verified. 
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CHAPTER II 



THE MATHEMATICAL MODEL 

A mathematical model will be designed to compute the 
probability of kill for an approaching aircraft penetrating 
an air defense complex composed of heavy antiaircraft gun 
batteries, light antiaircraft automatic weapon batteries, and 
surface to air missile batteries. The purpose of this model 
is to determine a relationship of altitude versus probability 
at each major terrain point through the air defense complex. 
The model is designed to be general. No specific weapon 
system parameters in existence are used in the calculations 
Values for the parameters are of the order of those found 
in actual systems. The general nature of the model was 
selected in order to permit its use with any system and to 
prevent the necessity of classification under security 
regulations . 

PROBLEM FORMULATION 

Since it is necessary to compute the probability of 
kill at each terrain point and at altitudes of every 100 
feet above terrain to 3 >000 feet, a concept for what is 
meant by a kill probability at a large number of specific 
points in space along a route of approach must be defined. 

The probability of kill by a weapon system is not defined at 
a point, probability is calculated from the number of 
rounds the weapon system is capable of firing while the 
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target is within range. The kill probability at any partic- 
ular terrain point and altitude, is defined to be the 

cumulative probability of kill resulting from the fire of all 
batteries of each type that bear on the target from each 
battery's maximum range to the terrain point in question. 

The amount of fire that each battery may produce is further 
constrained by terrain masking at the lower altitudes as it 
may occur. For example, at some terrain point k and alti- 
tude i, a line is drawn back along the route of approach 
to the beginning of the coordinate system or to the point 
of contact with a terrain formation, whichever is the shorter. 
Figure 1 is a diagram illustrating a two dimensional view 
of the terrain showing terrain points with altitude points 
above the terrain and the course lines associated with the 
(k,i)th point, representing the kth terrain point and ith 
altitude. In the figure, point F represents the 4th terrain 
point and the 1st altitude. Drawn from this point is a 
line back along the course line of aircraft flight as pre- 
viously described. In this case, the line intersects a 
terrain formation at point G. The model assumes that when 
terrain interrupts a continuous flight, the aircraft will fly 
over the obstruction and continue on at the prescribed 
altitude. The line from G to F will represent a portion of 
the aircraft flight path. If a battery were located at 
point L, a probability of kill will be calculated from the 
number of rounds fired if sufficient time permits the firing 
of any rounds from any battery that bears on a target flying 
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the course from G to F, Such a battery would be located at 
point L. The value of probability calculated here will be a 
number greater than or equal to zero for altitude 1 at 
terrain point 4. These numbers become part of some func- 
tional relationship existing at terrain point 4. Another 
example of discontinuous flight is at points H and I, The 
uninterrupted case is seen at point A, G, D, and E. Here 
no terrain obstruction enters into the situation and the 
factors effecting the probability calculated depend on the 
length of course lying in the detection radar beam pattern 
and the weapon system effectiveness envelope. As an example, 
consider a battery located at terrain point 5« If the 
model were evaluating the probability of kill contributed 
/ by this battery up to point (4,2), the number of rounds that 
could be fired as the aircraft traveled from point B to point 
(4,2) would be the determining factor, B is the point where 
the flight path passes into the radar beam pattern while 
point A is where the flight path enters the weapon system 
effectiveness envelope. The portion of the flight in 
both radar and battery effectiveness envelope only, enables 
firing to occur. Since the altitude of approach is so low, 
the aircraft cannot first be detected until point B. At 
point (4, 3) the battery at terrain point 5 detects the 
targets at point M In advance of its entry into the battery 
effectiveness envelope, A normal fire mission can be con- 
ducted in this case with adequate warning time. 
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Coordinate System 



In order for calculations to be carried out by the com- 
puter a coordinate system must be placed on the terrain 
section of interest. No existing map coordinate system 
will suffice. The coordinate square must be small enough to 
provide an acceptable approximation of the terrain contours. 
Figure 2 is a diagram from a top view of the terrain and 
battery situation. This view shows two dimensions, distance 
along the course line, and offset distance from the course 
line for terrain and battery positions. 

The grid system shown in Figure 2 shall be the size 
necessary to cover the section of terrain holding the air 
defense complex of interest. A route of approach to be 
analyzed is selected and the grid system is placed over this 
route on an appropriate map with the x axis parallel to the 
course direction and centered over the line. At the inter- 
section of each grid line, the altitude taken from the 
terrain contour lines is recorded along with the X and Y 
coordinates. The altitude above sea level becomes the H 
coordinate. All distances in this model are in feet. The 
X, Y, and H values are then punched on IBM cards for placing 
on computer tape as data input. The job of preparing map 
data for transfer to computer tape is tedious and time 
consuming. For an individual to attempt the task would be 
impractical. A team of trained operations plotters equiped 
with the proper equipment would be capable of accomplishing 
the task. Since no three dimensional data was available, 
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the computer model was solved using two dimensional data. 

The model was verified with this data except for the line 
of sight calculations for batteries not on the target flight 
line. The two dimensional terrain was taken from reference 
1. Each major terrain point was selected as either a peak 
or a valley. A sample of the terrain is shown in Figure 3 
and 4. Figure 3 is the detailed profile of the terrain. 
Figure 4 is the simplified version obtained by drawing 
straight lines between major terrain points. For a sample 
of the terrain actually used in the model see Figure 5 
through 7. The coordinates of the major peaks and valleys 
in the X direction and altitude in the H direction are re- 
corded in Appendix 2. The offset distance Y is zero for all 
computations but is carried along in the model for possible 
future use. Coordinates of the battery positions used are 
found in Appendix 2. Referring again to Figure 1, it should 
be noted that the X coordinate has its lowest value at the 
entrance of the air defense complex and increases to the 
maximum X coordinate value at the defended area. 

Air Defense Complex 

The three general types of weapons systems are the 
antiaircraft automatic weapon battery, the medium to heavy 
antiaircraft artillery battery, and the surface to air 
missile battery. These three are used since they cover the 
spectrum of ground to air weapon systems to be encountered 
in the reasonable future. 



19 



DETAILED T.ERRAI1I PROFIPR 




FIGURE 3 



SBtfLTFIFi) T3RRATU PROFIJJS 




FIGURE )(. 



20 



ALTITUDE ABOVE SEA ISVEL IN FEET 



l a 


O 

ir 


o 


O 


u 


. u 




21 



ALTI7HDK ABOVE SKA 1K7KL Hi FKB'i' 



vn (>. 

V •> u> 

8 8 




22 



ALTITUDE ABOVE SSA ILVi'.l III FEKT 




23 



The automatic weapons battery is defined to be a 10 gun 
unit of rapid firing weapons, five mounts, two guns to a 
mount. The fire control is by optics with target course 
and speed estimated by the crew. a computer is used to 
generate predicted azimuth and elevation for positioning 
the gun. Early warning is by radar with target acquisition 
by optical means. The medium to heavy antiaircraft artillery 
battery is composed of 4 guns with early warning radar, gun 
laying radar, and computer generated gun positioning signals 
via appropriate servo systems. 

The surface to air missile system is a low to medium 
altitude system. It is capable of detecting targets within 
its early radar beam pattern which is approximately from 1° 
elevation to 50° for volume coverage. The radar is assumed 
to possess effective moving target indication at the low 
altitudes. The tracking radars are also assumed to have 
low altitude capability. Since only one aircraft is attack- 
ing the complex in this model it is not necessary to specify 
the type of guidance with respect to the number of targets 
the system may attack at one time. 

The batteries are placed in the 2 dimensional terrain 
on the course line at various terrain points. The effect of 
terrain and increasing altitude, and the cumulative effect 
of more than one battery firing on a target is obtained from 
the various battery locations. No battery is permitted to 
fire at outbound targets. Each battery may fire from its 
maximum range or detection range up to its minimum range or 



24 



to its maximum turning rate in azimuth to the terrain point 
in question, whichever occurs first. No battery is permitted 
to fire on targets at a lower altitude than its own since 
this is not feasibly done. 

Neapon Target Kinematics 

In order to determine the length of target course line 
exposed to the fire of any battery, it is necessary to define 
the parameters of weapon target dynamics. Figure 8 is a 
three dimensional diagram of the problem to be solved. The 
line C D is the path the aircraft flies as it passes into 
the range of battery J as evaluated at the (k, i)th point. 
Line A B is the projection of target course on a horizontal 
plane at the same altitude as the jth battery. The .AH, the 
difference in altitude between the jth battery and the ith 
altitude of approach, is given by 



The difference in position of the kth terrain point and the 
jth battery along the course line is given by 





(1) 







( 2 ) 
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Likewise for the offset distance of the jth battery 







( 3 ) 



where 



A\; = V ' 00 . 




( 4 ) 



N = number of terrain points where a probability 



versus altitude function will be calculated 



30 = number of altitude points above terrain for 



which probability calculations will be accomplished 



= number of weapon batteries employed in 



the terrain. 

The other parameters are defined as follows: 

R = slant range from the battery to the 
s 

(k,i)th point 



R Xs = horizontal component of R g 

Rjjj = maximum effective range of the weapon system 
R ym = horizontal component of R^ 

R. = horizontal component of the range at which 



the lower portion of the radar beam, at the estimated half 
power point of the radar beam pattern, intersects the target 
course line 
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Xq = the distance J L along the course line from 
the maximum open fire range of the system to the (k, i)th 
point being evaluated. This parameter is the key distance 
which will determine the probability to be developed later. 

^DR = the warnin & distance, H J from maximum 
detection range to maximum open fire range 

© = the estimated elevation of the half power 

point of the lower edge of the radar beam pattern. 



Calculating Equations of the Model 

In this section all calculating equations will be 
explained or derived as appropriate. 

Prevent Outbound Firing . In a high threat operating con- 
dition when enemy attack is probable or imminent, air 
defense commanders are more concerned with the inbound raid 
than outbound raids. If a target passes through the fires 
of a particular battery, the probability of an outbound 
kill is less in most cases than an inbound raid. The danger 
of firing at outbound raids results from not being ready to 
engage the next inbound raid. Therefore the policy of not 
firing past the position of any battery is designed into the 
model. To insure no outbound firing £X is calculated in 
equation (2) without the absolute value notation. 



r*k/,-&*j > ax< ° 

O • A* >0 
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Since any terrain point with a greater X coordinate value 
lies in an outbound direction from the jth battery, AX > 0 
in this case. Firing is permitted up to the minimum range of 
the Jth battery, therefore setting AX = 0 permits later 
calculations for Xq up to B x j rather than to X k ^ . Any 
terrain whose X coordinate is less than B x j will produce a 
AX > 0 and bX will be set equal to /X ki - B x j/. 

Determine Slant Range . The slant range to the (k, i)tn 
point is given by 



R s - 



( 6 ) 



In Figure 8, it can be seen that R is the vector sum of 
AX, AX, and AH. Therefore equation (6) will hold. For any 
battery j whose maximum range R m < R g will not contribute any 

probability of kill. Therefore no calculations for the jth 

_ 

battery will be carried out in this case. 

Determine if Target Elevation is Negative . Since firing at 
elevations below 0° is not feasible for many reasons, this 
model excludes such cases. After solving for aH in (1), 
if AH<0, excludes the jth battery from a probability 
calculation. Figure 9 illustrates this case. 

Minimum Firing Range . In certain cases which arise in this 
model the (k, i)th point lies over the jth battery or is 
inside the dead zone of battery effectiveness. The model 
must determine when this situation occurs and permit firing 
up to the minimum firing range for the ith altitude of 
approach. Figure 10 is a diagram of the minimum range 
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Figure 9 



situation. Rp is the minimum firing range for the system. 
Angle B is the maximum elevation angle ivhich is the minimum 
of the maximum radar antenna elevation angle and maximum 
elevation angle of the gun or launcher. 

R X F = cos & ( 7 ) 



R as previously defined is given by 

xs 



( 8 ) 



from Figure 8. If R^ 3 < R^p, then firing must terminate at 
a horizontal range of R^p rather than B xs . Therefore equa- 
tion (6) becomes 



R 




^ F X + > 



> ^XF 
^ ^XF 



Maximum Azimuth Tracking Rates . When the jth battery is at 
some distance AY from the course line, it is possible for 
the maximum tracking rate of the system to be exceeded before 
the (K, i)th point is reached. If this should occur the 
distance Xq would be decreased. The amount of this decrease 
is calculated by finding that point on the course line where 
maximum tracking occurs and determining if it lies between 
(k, i)th point and the crossover point which is where aX = 0 
or before the (k, i)th point is reached. If the point of 
maximum tracking rate occurs before (k, i)th point on the 
course line then Xq is reduced by the distance between the 
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point of maximum tracking rate and (k, i)th point. Figure 
11 is a diagram of the azimuth tracking rate problem on the 
horizontal plane. Point K is the point where firing must 
terminate projected on the horizontal plane at the same 
altitude as battery J. The tracking rate must first be 
determined at this point. Let AZ be the azimuth tracking 
rate. Then 



A'Z 



= V XT 



# 



where V = tangential velocity of the aircraft at point K 

and R = horizontal component of the range to point K and 

xs 

was calculated in equation (8). From Figure 11 it is clear 



that V xT is given by 




VrAV 



Substituting this result into the equation above yields 



A'Z = 



V T XV 

i 



( 10 ) 



Now substituting equation (8) into (10) results in 



A'z - 



ay 



(ii) 



If AZy - maximum tracking rate for the system engaged, and 
AZ < AZ^, no constraint due to tracking rate is active. If 

however AZ > AZ T ., then the constraint is active and firing 

M 

must terminate at the point where the maximum rate occurred. 
To find this point solve equation (10) for R^ s when AZ is 
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( 12 ) 



replaced by A; Z M resulting in 



^5 




Using the value of R determined in (12) a new value for 
R s is found by solving equation (8) for R s as follows: 



Equation (13) has changed the (k, i)th point to reflect the 
constraint due to the maximum azimuth tracking rate. 

Maximum Elevation Tracking Rate . The same problem occurs 
in elevation as in azimuth for the possibility of exceeding 
tracking rates at low altitudes. In some cases the elevation 
tracking rate is exceeded at a range further from the (K, i )th 
point than minimum range. It is necessary to compute the 
elevation tracking rate and determine if it is an active 
constraint. Figure 12 is a three dimensional diagram 
illustrating the vector component of target velocity tangen- 
tial to the battery target line of sight. The tangential 
velocity in elevation is given by 




(13) 



V XT E = R 8 



(14) 



where the 




S 
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Since the velocity V^j is also opposite the angle A in the 
smaller triangle 



V XTE = V T sln A 



V 



XTC 

Using the relationship as before that EL = " — , results in 



HS 



EL 









( 15 ) 



If EL > ELj,j , where EL is the maximum tracking rate in eleva- 
tion of the particular system, then it is necessary to find 
the range at which the maximum rate occurs. Solving ( 1 6 ) 
for R ss gives 



R 



SE 



V 






( 16 ) 



The new value of Rg is determined from 
R s = MAX (R sa , R se ) 
and the new &X is determined by 

AX = V (1?) 

S 

At this point the conditions for terminating the firing 
at or before the (j, l)th point have been found. It is now 
necessary to find the point of open firing and then evaluate 
the effect of terrain on the length of course line between 
the open fire point, cease fire point, and the battery lo- 
cation. 
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Open Fire point. Referring to Figure 8, it is seen that the 



open fire point is at Ho The distance Xq is from the cease 
fire point which may or may not be the (k, i)th point to H. 
Depending upon the value of H, the distance X nD may have 
positive, negative, or zero values. Point f must be found 
as it is the point where the detection radar first intersects 
the target course line. Figure 13 diagrams the situation. 
Lines 1 and 2 represent different altitudes of approach and 
the different points of detection in range. The horizontal 
range to the point of detection is found as follows: 




4txn © 



( 18 ) 



This relationship will be a good approximation since the 
ranges where the beam pattern begins to curve are far beyond 
the maximum effective range of any weapon system. The hori- 
zontal component of the maximum range of the system is 
determined to be 

( 19 ) 

The distance X„ / X^n on the target course line must be 
D iJri 

above all terrain points and must be visible over its entire 
length from the battery position in question. The effect 
of terrain on low altitude approaches on the distance 
X^ / X may be such that its value is decreased to less 
than zero when no firing is possible. A method of evaluat- 
ing the terrain’s effect on the problem will be developed. 
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Terrain jgvaluation. The method to be used for determining 



if line of sight exists from the point of radar detection to 
the battery position is to compute the slope of the line 
from this point to the battery as well as the slope from 
each intervening terrain point to the point of radar detec- 
tion. If the slopes of the terrain points are greater than 
the slope to the point of radar detection then line of sight 
exists. If one or more terrain points extend above the line 
joining the battery and the point of radar detection, the 
line of sight does not exist for the length of the course 
line beyond the masking terrain points. A new point of 
radar detection is then found. A new slope line is com- 
puted from the masking terrain point at the altitude of 
approach and the battery. As before, all remaining inter- 
vening terrain point slopes are examined to determine if any 
more are masking the target course line. If another point 
has a slope less than the battery radar detection point 
line, then this higher terrain point becomes the new radar 
detection point. Depending on the terrain, the distance 
Xd / may be uneffected or reduced to zero. If X NK is 
the coordinate of the radar detection point then the slope 
to the battery position is given by 

c _ ~ 

^LPB 0^ -X NK (20) 

The slope of each terrain point to radar detection point is 
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(21) 



S 



LPm 



where 

Aj, . = altitude of (k, i)th terrain point 
h, 1 

H m ^ = altitude of the (m,l)th terrain point 
m = rath terrain point between the battery and 

the radar detection point 

1 = coordinate index of the offset distance 

from the course line 



^ = X coordinate of the intervening terrain 



points , 



As long as the engaged battery is located under the 

target course line, the determination of major terrain 

points that may mask the line of sight is a relatively easy 

matter and the slope calculations are simple. Figure 14 

is a diagram of the calculation of slope. Points 1, 2, 3> 

and 4 represent terrain points which may be typically found 

between battery j and the radar detection point L. As shown 

in the diagram, terrain 3 and 4 have greater slopes than 

the line B^,L. They do not mask the line of sight. Terrain 

point 2 has a smaller slope than E ,L and does mask the line 

Xj 

of sight. Therefore, point ^ is the new radar detection 
point at which a new comparison slope must be calculated. 

Let Z be the distance from ^ to B v - on the horizontal. Then 
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FIGURE lit 



( 22 ) 



Z - L ^ H l ~ ~~ X 2 j) 



This equation is derived from similar triangles in Figure 14. 
The value of the X coordinate of point 3 is given by 

= “xj _Z - , (23) 

When the engaged battery is located at some offset 

distance, Y from the target course line, a procedure must be 

found to determine if line of sight exists between the 

battery j and all the terrain that lies between the lines 

joining the cease fire point, and radar detection point on 

the target course line. Figure 15 illustrates this problem. 

Notice the coordinate grid. At each intersection of a grid 

coordinate, an IBM card recorded the X, Y grid index along 

with the altitude at the intersection. Notice also that 

the X coordinate is parallel to the target course line. 

As previously indicated, the X coordinate increases in the 

direction of target flight. Point J is the battery position. 

In order to determine if line of sight exists from J to X, TT , 

the slope of all terrain located on the line J, will be 

computed and compared to the slope of the line J, X MT ... In 

order to find the vertical and horizontal differential for 

all points on line J, X , a series of triangles will be 

NK 

solved for the hypotenuse such as a, b, c in Figure 15. The 
distance ac is the horizontal distance for the slope 
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FigurG_ 15> 



the X 



calculation. The distance be is, if added to X , 

NK 

coordinate of the closest point to p. Using the altitude 

of point p which has been recorded as previously described 

in digitalizing the terrain, the slope of the line X , C is 

NK 

calculated. The distance be is given by 

be = ab tan a (24) 



where 



ab = distance between Y grid lines 



tan a 





Y 



so that 



be- - 

Ay 

and 

<xc *■ \j (ab^ + ( bO 2 - 



(25) 



( 26 ) 



For each successive triangle, ab is replaced by Nab where 
N = (B - Y, ) / Yk - 1 

yj k »l 

and 



B = Y coordinate of battery i 
YJ 

Y = Y coordinate of the target course line. 
k,l 
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Substituting as appropriate, 



ol c = 



c ©ki y* r> 



7 K 






'Ay 



be* f Byj - Y K | N ) fSxj " ^Nk) 



( 27 ) 



(28) 



Solving (27) and (28) for each value of N, enables the slope 

to be found at each Y k interval and establishes line of sight 

along the particular line being examined. If any terrain 

slope is less than the comparison slope, then line of sight 

does not exist on that particular line from J to ^ . In 

this case X„ is increased by one X coordinate grid distance, 
NK 

X k and the slope calculation process is again repeated. 

The value of X which remains after each slope line is 
NK 

compared with terrain determines the actual radar detection 
point. In the three dimensional case it is possible to 
track a target up to a terrain formation masking a portion 
of the target course line between the cease fire and radar 
detection point, and either continue tracking through the 
mask in memory track, enabling relock, or to lose lock and 
require reacquisition on the other side of the mask. In 
the relock situation, the previous radar detection point 
X. ITr stands. If reacquisition is necessary, then a new 
value of X^v is computed from the point where the terrain 
ceased to mask the course line. The value that X. TT , takes on 
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is determined by 



T- f- 



T 



( 29 ) 



where 



T = time for aircraft to pass terrain mask 

X 



Xp = distance on course line masking line of sight 



to the battery 



V = aircraft speed 



If T v < T , where T is the memory tracking period, then 
-X- rn ni 

remains unchanged. If T v > T , X , receives a new value as 

X m’ nk 

shown below, 



mask 

Determination of Length of Course Line for Firing . W i t h 
the computation of X and the cease fire point after con 



sidering maximum effective range of the system, maximum 
radar detection range, and effects of possible terrain 
masking, the distance - X^ k is available for firing 

after the system delays have been accounted for. The system 
delays assumed in this model are as follows: 





( 30 ) 



inhere 



X k = distance between X grid lines 



and 



= number of grid lines spanned by the terrain 
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(1) Reaction time, t 



r 



which is the average of time 



for the radar operator to recognize target pips on the plan 
position indicator of the early warning radar set in the 
battery operations center, identification of the target and 
decision to engage the target. 

(2) Transfer action time, t , which is the time for the 

cl 

fire control operators to receive an assignment, lock the 
tracking radar on the target account for computer settling 
time, and gun crew reaction time to fire the first round. 

(3) Time of flight of the first round of the system to 
reach the target 0 

Average values for (1) and (2) are assumed but time of flight 
is calculated from system parameters. An equation for 
solution of time of flight for gun systems is given by 



If the k is known, (31) Is a useful formula. If a system 
P 

graph of time of flight versus slant range is available, a 
regression equation is also useful. This model uses a re- 
gression equation approximateing a typical system curve of 
slant range versus time of flight. The concept is to place 
the first round on the target at the maximum effective range 




(3D 



where 



R o £ = range to open fire 

V q = projectile muzzle velocity 

k = projectile drag 
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of the system. If sufficient time exists for the target 
to travel from the point of detection to open fire range, 
then the system is assumed to expend the time 

t = t / t / t 
T r r a f f 

prior to the targets arrival at the open fire range. If 
the detection range is less than maximum effective range, 
then the time t^, will reduce the time for firing. If on 
the other hand only a part of the time t^ is expended prior 
to open fire range, then some proportion of the t^, will 
reduce the firing time. 



D 



t 



'T T 



(32) 



where 

D. = the distance the target travels during the 
system delay time 

As a result of terrain evaluation, X , , as was seen, 

nk 

may have a new value. Therefore, the range to the maximum 
radar detection point is given by 



“'"A'/ 1 * 



(33) 



This relationship is illustrated in Figure 15. Since detec- 
tion range may have decreased to be less than R , the open 

AjU 

fire range, as a result of terrain restrictions, is now 

given by 



R xm = 



(3*0 



0 



4 ? 



Actually, R is the horizontal range component. The slant 

AJU 

range to open fire becomes 

r of = \1r* + ax.*- + Ay fc 

OF X xm , (35) 



The distance along the course line from the radar detection 
point to open fire is given by 




Now the distance on the course line involved in reaction time 
t T may be all within X QR if X QR > t^, or part of t T may be in 

X DR X DR < ^T° X DR ~ 0 then fc T wil1 reduce X D by the 

amount t^V^o X^ will also be reduced by the time of flight 
for the last round fired since any round fired after the 
point (k, i) is reached does not contribute to the proba- 
bility of kill up to point (k,i)„ Therefore, Xq is given by 



^ - AX-VV'V; 

X Q ^ CKXp,<^37) 

- Ax-+ f \^r ; X PR ^ x T 

where 



X T1 X T ’ X DR 

x T = (t r + t a ) v T 
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Having derived the equation for X Q , the length of course 
line available for firing, it only remains to determine the 
number of rounds that will be fired in the time the target 
passes over X Q . The time for target passage is 



T = 






V- 



(38) 



Then the number of rounds fired is given by 



N 



r 




(39) 



where t . = time to fire one round and the rate of fire 

1 Vr 
is t 1 - • 

Vr 

If < X^, then the entire delay time including the time 
of flight of the first round occurs during the early warning 
time period and (39) becomes 

1- 4- l 

M r ' '/ r 




) 






( 1 » 0 ) 



> X 0R ^ X T 



Probability Model . In Chapter I, it was stated that the 
model was being developed for three different air defense 
systems. In the development of the model, the general 
scheme was to permit all calculations to apply to each type 
of weapon system except for the parameters used in the com- 
puter model. Except for the regression equations used for 
determining time of flight, this generality has been success- 
ful. However, the models must be separated for probability 
calculations into a separate one applying to each weapon 
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system. In the case of the two types of gun systems, a 
common model is still possible with one or two exceptions. 
The radar controlled antiaircraft gun system will be delt 
with first. The range at which each round is fired is 
computed as follows: 



This model assumes that the area a gun fires into is defined 
by a circular normal probability curve. The square of the 
radius of the circle described thusly has been determined 
empirically in references 2, 3, and 4 is given by 



<r± is the result of all the errors introduced by the fire 
control, data transmission, orientation, and alignment of 
the battery. Equation (42) is based upon the following 
theoretical equation: 



R,; =\|cx 0+ Ax-v T t /r '0 > + ^AV 3 - 



(41) 



1 = ‘ > N r 




(42) 



where 



= angular error of the gun in milliradians. 
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where 



= quasi static errors consisting of 

(1) Boresighting and alignment 

(2) Range errors (radar developed), radar 
timing errors due to mechanical alignment 

(3) Servo lag errors 

(4) Computer error 

(a) static errors 

(b) ballistic match errors 

(c) prediction errors 

(5) Battery emplacement errors 

S t = Tracking perturbation 

S„ = Plight roughness 
f r 

The constants and results of using equation (42) have been 
verified in actual test firings according to reference 2. 
Target Area , The average cross sectional area of the tar- 
get’s vulnerable sections is denoted by A in square feet. 

The vulnerable cross sectional area changes with different 
positional attitudes. An average value is assumed in the 
following calculations. In addition a method for calculating 
the area as a function of target attitude was also imple- 
mented. This area is computed from 







(43) 
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where 



A. = percentage of area presented for the ith 



round fired 



= the range of the ith round 
X. = the X component of the range to the target 



at the ith round fired 



A^. = cross sectional area of the standard target 
defined to be 400 ft^ in the X coordinate. 



800 ft' 



150 ft 



A = same as above for the Y coordinate set at 



A = same as above for the H coordinate set at 
H 



Y^ = the Y component of the range to the target 



at the ith round fired 



H = the H coordinate of the range to the target 



at the ith round fired 



The standard target may then be compared with any other 
target for converting probabilities. 

The probability that one round from the gun hits the target 
is the ratio of the vulnerable cross sectional area of the 
aircraft to the area the weapon fires into as follows: 




A 



A7Tcr: 



> 



) 



♦ 






♦ «• « ^ N R. 



(44) 



The probability required is the probability that one round 
kills the target. If the probability of kill given a hit 
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is P[K/Hj then the probability of kill, assuming indepen- 
dence, is P[K/Hj p[h]. Using this fact equation (44) 
becomes 



P. 



A 

^ rr<r 



PC K 



(^ 5 ) 



If N_^ is the number of rounds fired, assuming independence 
of rounds fired from the gun, the probability of kill is 
given by 




where 




(46) 



N - number of guns in the battery, 

S 

guns firing the probability that the aircraft 
assuming independence of guns in the battery, 

P K ; - '-(l-PO^ 



With all 
is killed, 
is given by 



(4?) 



The second type of weapon system is the automatic weapon 
battery. This battery does not possess gun laying radar. 

It depends on the optical means for acquiring the target. 
Early warning radar is assumed available. All equations in 
the previous model are applicable except equation (42), 

The total error for a gun with optical sights is assumed 
circular normal with the square of the radius of the circle 
of error given by 

*T “ (48) 

» 
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This equation like (42) has been empirically determined and 
verified according to reference 2, 



In the surface to air missile battery, the third type 
of weapon system, a simplified model is assumed. Specific 
models are not feasible in order for the computer model to 
be applicable to any type of missile system. Therefore, the 
single shot kill probability is p. Then the battery kill 
probability is given by 



N was calculated previously. However, in the missile 



battery the number of rounds fired at any target is con- 
strained by the basic load of missiles. Further the high 
single shot kill probability of a missile round precludes 
the firing of more than 2 or 3 missiles per target in most 
cases. This model insures that the value of N permitted 
is as follows: 



Any particular battery firing policy may change the value 



culated for all batteries. Then the overall kill probability 
for all batteries is given by 




(49) 



r 






r 



# 



(50) 



of as the tactical situation may dictate. 

For each terrain point k and altitude i, P is cal- 

kj 
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i — I j A , * ♦ » > 30 

j = N, 

4 

The assumption of independence is quite valid for equation 
(51) since all batteries act independently in actual firing. 
Equation (46) also requires the assumption of independence 
between each round fired. If the weapon is a radar con- 
trolled gun using equation (42) to model the area within 
which all rounds are fired, the probability of a hit for 
each round , must be independent for each value of i for 
the assumption to hold in (46), When examining the factors 
which determined (46), it is seen that certain of the errors 
are the same from round to round or are corrolated errors. 
These errors tend to break down the assumption of inde- 
pendence. Most of the errors are random in nature from 
round to round and are the factors supporting the assumption 
of independence. It is not known how serious the correlated 
errors effect the results of equation (46). If independence 
were not assumed, (46) would be modified as follows: 
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The difference between the probabilities in (46) and the 
conditional probabilities in ( 52 ) constitutes the error if 
if is significant. All models found in the literature use 
the independence assumption. Investigation of the proba- 
bilities in ( 52 ) to determine how bad the assumption of 
independence really is may prove interesting. 
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CHAPTER III 



THE COMPUTER MODEL 

In order to perform the calculations necessary to 
produce a relationship of altitude versus probability at 
each terrain point within the range of any battery in the 
defensive complex, a computer program capable of performing 
the calculations is essential. If the resulting computer 
program is to have any lasting value, the program itself 
must be documented in order to make clear how the program 
was designed and how the problem was solved by the computer. 
An explanation of the Fortran variables is found in Appendix 
III with the Fortran program located in Appendix IV. 

The program was originally designed to produce as an 
output an array of probabilities and corresponding altitudes 
for each terrain point examined. It was expected that the 
results would be punched on an output deck for use as an 
input to the computer program in reference 1. The output 
deck from this program would then replace the probability 
calculations in the other program. This replacement has 
not been possible for reasons to be explained in a later 
chapter. 

Computer Inputs 

The inputs to the computer are of two types. The 
terrain digitalization and the weapon system parameters. 

The three arrays X , Y , and H are the coordinates 



5 ? 



of the major terrain points being examined. In the two 

dimensional problem this is all the data needed. If the 

three dimensional data is available, three more arrays 

would be needed to record the terrain contours through out 

the area of interest. The arrays HH for altitude, XX, , 

I , K I ,K* 

X coordinate, and YY the Y coordinate of each recorded 

I»K 

terrain point would also have to be inputs to the program. 

As previously explained, the model was developed for 
any weapon system. If the model is to represent any specific 
weapon system the input parameters necessary to define the 
system must be clearly understood. An explanation of the 
input parameters is contained later in this chapter. The 
parameters are defined in the Fortran code notation to 
facilitate their incorporation into the computer program. 

Description of Computer Solution . 

Figures 1 6 , l6a, and 17 are a block diagram of the com- 
puter program in general terms. The following description 
will relate to the block diagram. 

Block #1 . All data is read into the computer with units in 
feet and seconds ,, angles are in degrees. 

Block #2 . Most parameters of angular measure are given in 
degrees. The program converts all degrees to radiams. 

Block #3 » The program is solved once for each batteries 
at each of 30 altitudes and for each of N terrain points. 

This constitutes 30 I^N iterations which in the sample used 
in the actual program is 30 x 11 x 53 = 18,490 iterations. 
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Block Diagram of Computer Program 



h 




Figure 16 
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Block Diagram Continued 



12 



13 




Figure 16a 
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Block Diagram Subroutine LOS 
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Block #4. In order to permit the program to use the same 



code for solving the problem for all three types of weapon 
systems, the system parameters are converted to a common set 
which represents the values corresponding to the particular 
system being examined at any one time. 

Blocks #5, 6, 7 i 9 . Antiaircraft fire has been con- 
strained by 

(1) Firing only at inbound targets 

(2) Not firing at targets out of range 

(3) Not firing at targets lower in elevation 
than the battery 

(4) Not firing at less than minimun range 

(5) Not exceeding maximum azimuth and elevation 
tracking rates. 

A section of computer code as shown is Figure 16 provides the 
calculations to determine, at each iteration, the occurance 
of any one of the five constraints listed above. The overall 
effect of each of the constraints will now be explained. 

(1) The program is designed to permit firing up to 
minimum range or to the point of maximum azimuth and ele- 
vation tracking rates. If a terrain point is well behind 
the battery or in the opposite end of the course from the 
direction of target attack, the minimum point is thusly 
established. On the other hand, if the terrain point in 
question is in the direction of attack, then the terrain 
point itself is the cease fire point unless the above 
constraints act to increase the cease fire range. 
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(2) Firing at targets that will always be out of range 
is prevented by assigning zero to the probability of kill if 
the range to the terrain point and altitude is beyond system 
range, 

(3) If a target is lower in elevation than the battery, 
zero is assigned the probability of kill for that battery 

at that altitude, 

(4) Jhen either minimum range, maximum azimuth and 
elevation tracking rates are exceeded, whichever occurs 
first, the cease fire range is increased to that point. 

Block #10 , The subroutine LOS is called at this point which 
calculates the distance of course line available for firing. 
The function of LOS will be explained later. 

Block #11 , The length of target course line available for 
firing resulting from terrain considerations may produce a 
new value for open fire range. Based upon the value of X^, 
the open fire range is calculated. 

Block #12 , Depending on what type of battery is engaged, 
the program will go through one of the three paths in this 
block. The probability of kill depending on the number of 
rounds fired is calculated for one of the three systems. 

Block #14 , The overall kill probability is computed here 
for one iteration which includes all weapons at one terrain 
point and one altitude. 

Block #15 ° After all iterations are completed the program 
is capable of plotting the curves of altitude versus 
probability. Any other output is available when desired. 
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Subroutine LOS 



Block #1 . Receiving all necessary data from the main pro- 
gram, all calculations planned are possible. Radar warning 
time is determined prior to terrain evaluation. The value 
for the point of first detection is the result of this 
block. 

Block #2 . Here the program determines which terrain points 
lie between the battery and the point of first detection. 

If any terrain interrupts the line of sight, a new value 
for radar detection point is calculated. The computations 
are for two dimensional terrain which mean terrain for 
batteries on the target line of sight. 

Block #3» If three dimensional terrain is used, sub- 

routine SEE is called where the point of radar detection, as 
a result of terrain, is found. 

Block #5 « Open fire range without regard to delay is 
calculated for use in determining time of flight. 

Block 6 . Again the program branches into three paths, one 
for each weapon system. This is necessary since each system 
has a different function of slant range versus time of flight. 
Block #7 . The delay time which is the sum of crew reaction 
time, system delay time, and time of flight of the first 
round fired is calculated in this section. This value is 
not used in the case where early warning time is less than 
total delay time. 

Block #8 . After all the above calculations are completed, 
it is possible to find the actual point of open fire. This 
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point is the range at which the first round is placed during 
delay time to maximize the number of rounds fired. Three 
possible situations arise in calculations for the open fire 
point as shown in Chapter II » equation ( 35 )« These are 

(1) Adequate warning time exists for all delay time 
to occur during the early warning period, 

(2) Only part of the delay time occurs during the early 
warning period and the open fire point is reduced, 

( 3 ) All of the delay time must fall during possible 
firing time due to the radar detection point being less 
than maximum effective range. The subroutine sends 

back to the main program . 

A Description of Input Parameters 

At some later time, it may be desirable to compute 
the results from the model for an actual air defense system. 
In order to do this, an explanation of the input parameters 
is essential. A description of how to insert new inputs 
is provided below. Statement number 4 in Appendix III is 
the data statement. Contained in this statement are all 
the system parameters. The program is convertable to any 
specific weapon systems by changing the values in the pro- 
gram data statement, to conform to those of the new systems. 
Included with the description of each parameter is the value 
assumed in the current form of the program. 
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Parameter 


Description 


NAP 


Number of altitude points 
above terrain examined. 
NAP = 30 


NTT 


Number of terrain points 
examined, NAP = 84 


NFU 


V 

Number of batteries placed 
in the terrain. NFU = 11 


ELRMG 


Maximum elevation tracking 
rate for the radar controlled 
gun system. ELRMG = 10°/sec 


ELRMA 


Maximum elevation tracking 
rate for the automatic weapon 
battery. ELRMA = 20°/sec 


ELRMM 


Maximum elevation tracking 
rate for the missile battery. 
ELRMM = 15°/ sec 


AZRMG 


Maximum Azimuth tracking 
rate for radar controlled gun. 
AZRMG = 12°/ sec 


AZRMA 


Maximum azimuth tracking 
rate for automatic weapon 
battery. AZRMA = 20°/sec 


AZRMM 


Maximum azimuth tracking 
rate for the missile system. 
AZRMM = 15°/sec 


RMG 


Maximum effective range for 
radar controlled gun battery. 
RMG = 45,000 feet 


RMA 


Maximum effective range AAV/ 
battery, RMA = 6,000 feet 

1 V) 


RMM 


Maximum effective range 
missile system. 

RMM = 105*600 feet 

-is,l 


THETI 


Elevation angle for lower 
portion of detection radar 
beam pattern at half power- 
point. THETI = 1.1°2 
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Parameter 


Description 


THET 2 


Elevation angle for detection 
radar, A AM battery, THET 2 = .5 


THET 3 


Elevation angle for detection 
radar, missile battery, 

THET 3 = 1,°0 


API 


The constant term in the linar 
equation for slant range versus 
time of slight, gun battery, 

API = 0 


CONT 


The constant in the expo- 
nential equation for time of 
flight for AAM Battery, 

CONT = ,0002 


YK 


Distance between Y grid lines, 
three dimensional terrain, 

YK «= 500 feet. 


VT 


Target speed, VT = 350 ft/sec 


TRG 

TRA 

TRM 


Crew reaction time. This time 
is measured from time of first 
decision of detection of tar- 
get to track radar lock on. 
This is an average figure. 

TRG = gun system, TRA = AAW 
system, and TRM = missile 
system, TRG = 30 sec, 

TRA = 10 sec. TRM = 30 sec. 


ATG 

ATA 

ATM 


System delay time measured 
from radar lock on to first 
round fired for gun battery, 
AAM battery, and missile 
battery, ATG = 20 sec. 

ATA = 5 sec. ATM = 10 sec. 


LT 1 


Total number of gun batteries. 
LT 1 = 4 


LT 2 


Total number of AAV/ batteries 
plus gun batteries. LT 2 = 10 


LT 3 


Total number of gun batteries 
plus AAN batteries plus 
missile batteries. LT 3 = 11 
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A 


Cross sectional area of target. 
The average value as target 
aspect changes. A = 350 feet 


SIG 


Angular error of radar^ con- 
trolled -AA gun when firing at 
a fixed point. SIG = 1.5 raille 
Tadians 


GUNS 


The number of guns per gun 
battery.' GUNS = 4 


TCRA 


Cyclic rate of fire, AAW 
battery. TCRA = 6 rounds/sec 


PKSM 


Single shot kill probability 
for missile system. PKSM = .48 


NR 


The firing policy of the 
missile battery, NR rounds 
per target. NR = 2 


TIM 


Time for one round fired from 
the gun battery. TIM = 20 sec 


AGUN 


Number of guns per AA«^ battery, 
AGUN = 4 ' 


TX 


Memory track period for all 
systems. TX = 15 sec 


LY 


Index number for the Y coor- 
dinate of the target course 
line terrain points. LY = 1 


RMPG 


Minimum firing range,' gun 
battery. RMPG = 500 feet 


RMPA 


Minimum firing range for AAA 
battery. RMPA = 500 feet 


RMFM 


Minimum firing range for 
missile system. RMPM = 9500 
feet 


PE 


Maximum elevation angle for 
all track radars. PE = 87° 


YM 


Slope of the linear regression 
equation for time of flight of 
the gun battery. YM = .0007 
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ZM 



Constant in linear regression 
equation for time of flight 
of missile system. AM = 2,000 



SM 



Slope of linear regression 
equation. SM = 1200 



PG 



Probability of a kill given a 
hit for the gun battery. This 
Is the lithality factor. 

PG = .8 



PA 



Probability of kill given a 
hit for the AAN battery. 

PA = .3 



Certain other parameters may be found in the data list. If 
they are not included in the list above, they are not in use 
in the final form of the program and should be disregarded. 

Terrain Data Input . Statement number 7 in the main program 
is the read statement which picks up the data punched on 
84 data cards. Here the variable NTT is set equal to 84. 

If in a subsequent use of the program more terrain points 
were required, changing the value of NTT and providing the 
same number of data cards provides the change. Statement 
number 8 reads in the coordinates of the Batteries. It is 
essential that the first system type be radar controlled 
gun batteries up to LT 1. Then the next type must be AAri 
batteries up to LT 2 which was defined to be LT 1 / the 
number of AAW batteries. The last system in order is the 
missile battery. Generally only one missile battery will 
be in range of a single azimuth of approach unless the 
azimuth is in the overlapping zones of fire of two adjacent 
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batteries. Setting the values of LT 1, Lt 2, and LT 3 along 
with the proper ordering of the data cards of each type of 
weapon system constitutes the conversion to a new defensive 
complex. 

Three Dimensional Data Input . If three dimensional data is 
to be used, three new arrays must be introduced for the X, 

Y, and H coordinates of each recorded terrain point. //hen 
this is done, statement 6 in subroutine SEE must be re- 
placed by a statement using the new altitude array for 
H(J,ly) presently in the program. 



70 



CHAPTER IV 



PROGRAM VERIFICATION, RESULTS, AND ANALYSIS 

The underlying factor in the model design from concep- 
tion has been that the air defense capability is maximized. 
This means that the attacking aircraft situation is the 
worst possible condition for survival. This fact is appar- 
ent when considering the general model structure presented 
in Chapter II. It should be recalled that the aircraft 
was required to penetrate the defense from out of range to 
the defended area flying a straight and level course, with 
a raid size of one aircraft, and no electronic counter 
measures. In addition, the defense was assumed to be alerted, 
manned, and concentrating on the single intruder. The air- 
craft survivability index calculated in the model when 
submitted to a programing algorithm for minimization within 
the aircraft flight constraints finds the flight path' which 
minimized the probability of kill under the worst possible 
attacking conditions. In actual application, anything the 
aircraft may do relieving the sever restrictions, wholly or 
in part, increases the survivability index. The act of 
following the programed output recommended flight path, 
penetrating the defense with more than one aircraft in the 
raid, and using • ele-ctronic counter measures will increase 
survivability. These facts must be considered when evalu- 
ating the model’s results. 
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tfith the development of the model and the program for 
the computer, a verification of its main computing aspects 
is in order. Following this, the computer results will be 
presented and analyzed. 



Verification of Computer Results 

The two main problems solved by the computer within 
the mathmetical model are the calculations for determining 
line of sight, in two dimensional terrain only, and the 
calculations for the probability of kill given a positive 
course line. All other computations within the program 
support these main requirements. 

Terrain Evaluation . In order to verify the computer solu- 
tion of the line of sight, a sample of computer print out 
contained in Figures 18 and 19 is presented. Figure 18 is 
a case where the coordinates of the point of maximum radar 
detection, X ^ = 26l40. The range of first detection is 
163,060 feet since the coordinates of Battery 6 are given 
as B x £ = 189,200 on the X axis and 67 OO feet altitude*. The 
comparison slope S^pg = . 00873 » kk = 46, and II = 5« II is 
the terrain point number just greater than X n ^» KK is the 
terrain point number just less than B x ^. Thus, the program 
has determined the number of intervening terrain points to 
examine for slope comparison. As explained in Chapter II, 

if any slope S is less than then line of sight as 

LP LPB 

corrected by the distance z exists. Observing Figure 18 

again, it is seen that S (5) through S (29) are greater 

LP LP 
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Figure 19 
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it is found 



than s L pg» When S^(30) is compared to S^g, 
to be less and thus masking line of sight. Z = 7 0338.14 
was calculated from equation (22) Chapter II. A new value 
for X nk is given as X nk = 118862 from equation (23). Using 
the assumption that targets fly over the terrain when the 
flight path passes through terrain, a new comparison slope 
is calculated and X ^ is changed from 26l40 to 118862. 

Three more terrain obstructions exist such that the final 
value of x nk = 183740. From equation (37), Xg =-1953.7 
which indicates that no time is available to fire any 
rounds and battery 6 contributes no probability of kill at 
terrain point 62, and altitude 7. Figure 19 illustrates a 
case where terrain masking occurs once but a distance of 
4830.9 feet for Xg remains for firing remembering that the 
time for firing the first and last round have already been 
removed from Xg. Therefore three rounds are fired, one 
more than shown on Figure 19 which is the last round fired. 
The slope comparison calculation is performed 18,490 times 
less those cases where negative elevation angles occur and 
the terrain point is out of range of the battery being 
evaluated. 

Probability Calculations . A sample of computer output of 

probability calculations is found in Figures 19» 20, and 21. 

Figure 19 is a sample of a radar controlled gun battery. 

As was seen, three rounds were fired with a probability of 

kill of P (1) = .46. This is the battery kill probability 
k 

resulting from 12 rounds fired, 3 P er gun at a range of 
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Computer Probability Calculations 
Automatic Weapon Battery 
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Computer Probability Calculations 
i Gun Battexy 
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5089, 3059, and 1280 feet. Due to the short range, the p 
are .04, .08, and .1. Looking at Figure 21, a run of 17 
rounds was fired. Here the altitude eliminates terrain 



masking with a course length of X^ = 33762.9 feet. However, 



in the earlier example. The short range in the other example 
accounts for the high probability for only 12 rounds fired. 
Figure 20 is the calculations for an automatic weapon battery. 
Fifty rounds were fired in this case. Actually more rounds 
could have been fired but the constraint on this weapon 
system limits firing to 50 rounds. The range is about 4500 
feet on the average with p varying from .000016 to .0000217. 
The battery kill probability is p^(5) = .0064. 

Program Results . Appendix V is a set of sample curves which 
plot altitude versus cumulative probability of kill for 
specific terrain points. The program prints out the curve 
for each terrain point examined. Included in the set are 
illustrative examples of program results. 

Appendix VI is a set of curves plotting probability 
cost index versus altitude above terrain. This curve is 
obtained by the following relationship: 



the range has increased to 23,526 feet so that the P^ vary 

from .00046 to .06l. The battery kill probability, 

P (4) = .46. This probability is the same as for battery 1 
k 




(53) 



k - \ ) n « ^ Nl 

I : l N ' * ■> 30 
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where 



P^ = the probability that the aircraft pilot 
detects the navigational target. 

The values for P^ were obtained from the computer program 
VISTRAC in reference 1. 

Analysis of Results 

The purpose of the mathmetical and computer models of 
the air defense complex is to compute and plot a curve of 
provability of Kill versus altitude at each terrain point. 
Having obtained these curves, it is necessary to analyze 
them in order to determine the reasons for their shape. 

With this information, the application for the model's 
output can then be considered. 

Terrain Effects . One of the most significant contributors 
to the shape of these curves is the terrain within which 
the air defense is situated. As was seen in Chapter II, 
Figures 5 through 7> the terrain has different formations 
producing different curves. The curves in Appendix V 
illustrate this fact. 

Due to inaccessibility, many prominent terrain points 
are not occupied by firing batteries and only cause serious 
masking to batteries located on accessible terrain with 
less altitude than some others. - The result is that at low 
altitudes of approach, some batteries may be masked and 
contribute no probability of kill to the target. Whereas 
at higher altitudes the target is unmasked soon enough to 
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be fired upon by the battery,, Figure 22 illustrates this 
concept. The line of sight to altitude 1 intersects the 
course line at A and X^ = BE. From this result it 
is clear that the probability of kill at altitude 2, 

P Kk2 > P Kkl since *d 2 > X D1 due to incre asing the altitude 
of approach. With the effect of terrain masking above, 
more rounds will be fired as altitude is increased. 

It can be observed in Appendix V that some curves 
appear to have discrete jumps or regions of discontinuity. 
This is caused by terrain effects as well. These jumps 
occur, when at the loxNrer altitudes, a particular battery 
with high single shot kill probability is masked. At the 
mext altitude of approach considered, the battery is sudden- 
ly unmasked and a large jump in probability is the result. 
Effect of System Parameters . The battery effectiveness 
envelope is another influencing factor to curve shaping. 

This envelope was shown in Figure 1, Chapter II. In Figure 
1, notice the curve of the envelope as altitude increases. 
This curve is an assumption which approximates the effect- 
iveness of the average type system in use. If this model 
is used to simulate a weapon system with a radically differ- 
ent envelope, then a modification may be in order. The 
effect of the assumed envelope is also illustrated in 
Figure 22. At altitude 3> the line of sight intersects 
the target course line at G, but is beyond the range of the 
gun system. Therefore, = CF since C is the point where 

the target course line intersects the weapon system envelope. 
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Figure 22 



It is obvious that CP < BP by the amount of BC* from which 



it is determined that X^ < X^ and accordingly < ^Kk2' 

This then illustrates why the curves tend to slope down as 
altitude is increased over some regions of altitude. 

Another system parameter which effects the output 
curves is radar performance. In Chapter II, the half power 
point of a radar beam pattern is assumed to be the con- 
servative limit of usable volume coverage provided. There- 
fore, any target in this model falling within the beam 
pattern is assumed detected with probability 1. This 
assumes a lateral range curve to be 



P(b 



- £ 



(54) 



where 



FM 



r*> ) * 

o ; 



(55) 



Substituting (55) into (54) results in 




I \ X — R- 

0 ; x 



(56) 



These equations are found in reference 11. Figure 23 is a 
plot of such a lateral range curve. The model assumed no 
maximum value for R as a parameter input. The maximum value 
R can take on is limited by the length of the terrain coor- 
dinate system. Maximum R used for the computer program is 
47 miles. A nominal value for battery acquisition radars is 
65 miles. If a coordinate system were longer than the 
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Lateral Range Curve 
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Figure. 23 



83 



maximum rated range of a particular system, then a value of 
R for radar detection would have to be assigned. further, 
if the radar detection range normally used in the program 
were at the maximum rated range, the assumption of equation 
(56) may introduce significant error and a different radar 
detection rate would have to be incorporated into the model. 
The reason this factor was not developed was due to the 
extremely low altitudes of target approach where a detection 
range does not become greater than a value much less than 
the maximum rated range of the radar. Such a development 
would constitute a needless complication. 

Two other parameters x^ere used to perform a sensitivity 
analysis to determine the effect of a change in these 
parameters on the probability curves. The minimum and 
maximum range of the AAW battery was originally 500 feet 
and 6000 feet respectively. These values were changed to 
1500 and 9000 feet for one run of the computer. Curve 
number 8 in Appendix V is the result of using both para- 
meters plotted on one graph. Curve I is the result of the 
increased parameters. These two curves are not radically 
different, but the increased range did cause a general 
increase in probability for most values of altitude. The 
rate of fire and total number of rounds fired per run per 
gun is a more powerful factor in the model. Curve 5 in 
Appendix V is a curve for AAW batteries firing up to 200 
rounds per gun per run at 1000 rounds per minute with 10 
guns firing independently per battery. The result of this 
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situation is a very high probability curve which is obvious 
since if sufficient rounds are fired high probabilities will 
result. Practical constraints will not permit this type of 
performance in the usual case. 

Battery ^placement Effects . The probability curves gener- 
ated by the terrain in the computer model are found to be 
quite different from terrain point to terrain point. A few 
curves from adjacent terrain points are similar in the case 
where terrain does not change rapidly. Except for this case, 
the curves have different shapes. A different curve will 
obviously result if a different battery emplacement scheme 
is implemented since new battery to terrain masking situa- 
tions occur. As the program calculates probabilities for the 
different terrain points, new situations occur such as a 
battery on a high peak is in range or a battery behind a 
peak is unmasked. This variation in terrain conditions is 
why the curves are not exactly related even though the 
probabilities calculated are cumulative. The curves for 
terrain points 47 and 69 , curves 3 and 4, Appendix V, when 
compared, are found to be quite different. An additional 
factor to be kept in mind is that the altitudes above 
terrain are not at the same absolute altitude. Altitude 
600 at terrain point 47 is not necessarily the same altitude 
at terrain point 69 . There is no intention that there be 
any relationship between the different terrain points. The 
probability curve for each terrain point is separate and 
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distinct from any other unless, of course, ther is little 
terrain change between them except the X coordinate distance. 

Model Application 

In Chapter I, it was pointed out that the work done 
in reference 1 found the relationship between probability of 
radar detection and altitude above terrain to be approxi- 
mately linear. As a result, the computer program incorpor- 
ated a linear fit to the data. Following this the proba- 
bility of pilot navigation was added to the linear function. 
The functional relationship was necessary for use in the 
simplex algorithm for linear program solution for the optimal 
altitude of approach. At the beginning of this work, the 
concept was to develop the air defense model as an input 
to the program in reference 1. After examining the output 
curves of Appendix V a linear fit did not appear feasible 
for the data. In order to determine how bad a linearity 
assumption would be, a regression equation was calculated 
from the data of terrain point 62 in Appendix V. Curve 
number 4 is a plot of the equation. The function derived 
is given by 

Y = .128 + .000369 X (57) 

To submit the output data to a linear fit would Introduce 
serious misrepresentation of the model capability. This 
fact is obvious from the results of the linear fit to the 
data of terrain point 62 . 
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Probability Cost Index . The curves in Appendix VI are a 
probability cost index. The values are not true proba- 
bilities since they are the sum of two probabilities as 
shown in equation ( 53 )* The term cost is defined to be a 
measure of the material effort necessary to penetrate the 
defense, reach the target, and deliver the ordnance. Since 
either one of the events constitutes failure in the mission, 
the higher the value of the cost index the greater the 
chance that the mission will be a failure. Consequently, 
the survivability of the aircraft in the penetration effort 
is maximized when the cost index is minimized. 

ihen examining the curves in Appendix VI, it is ob- 
served that each curve contains a global minimum. The 
curves then possess a point which minimizes the cost index 
and provides the optimal survivability index number at the 
associated altitude. 



87 



CHAPTER V 



CONCLUSIONS AND EXTENSIONS 

The development of the air defense model within the 
assumptions and limitations is complete. As a result of 
the output from the computer model it is concluded: 

(1) That the probability calculations developed a 
result that varies with terrain, air defense systems, and 
employment tactics. Consequently, no assumption of linear- 
ity between probability of kill and altitude above terrain 
is justified. In fact, an investigation of the degree of 
polynomial that a least square fit would require to model 
the curves in Appendix VI is necessary. Such an investi- 
gation would be required to discover how the degree would 
vary with terrain points using the same system parameters. 
It would also have to determine how the degree would vary 
with different air defense systems and terrain. Before 

an algorithm can be found to solve for the minimum point 
on the cost curve, the curve itself must be functionally 
modeled. 

(2) That no practical results would be obtained from 
inserting the output from the air defense model into the 
program in reference 1, If such an attempt were made, the 
simplex algorithm used in reference 1 would obtain a solu- 
tion, but it would not be known if the solution were the 
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global minimum. Further, the linear regression curve would, 
grossly missrepresent the probabilities calculated, by the 
air defense model, 

( 3 ) That further investigation is necessary to discover 
the feasibility of fitting the curves in Appendix VI to 

some degree polynomial such that the degree chosen is 
acceptable for the curves of all the terrain points. If a 
single degree polynomial is not feasible, then a method of 
finding the best fit for each curve over a specified range 
of altitudes must be adopted. 

( 4 ) That having determined the functional relation- 
ship for the probability cost curve, a non linear programing 
algorithm must be found that solves for the global minimum 
of each curve within the aircraft flight path constraints 
used in reference 1. 

Air Defense Model Extensions 

Two major assumptions in the air defense model should 
be investigated further. The assumption of independence 
between rounds fired successively from the same gun may 
introduce serious error. Research into the determination of 
the relative weights given to correlated and uncorrelated 
error should give an indication of how unrealistic the 
assumption of independence really is. If independence is 
unacceptable then means for finding the conditional proba- 
bilities associated with dependently fired rounds should 
be studied. 
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An additional problem for further study is the assump- 
tion that the gun fire shot pattern is circular normal. 

Much evidence points to the fact that the shot pattern is 
actually elliptical normal. Empirical equations for the 
elliptical normal error pattern such as equation (42) and 
(48), Chapter II in the circular normal case should be used 
if such equations exist. 

The final extension is in the area of three dimensional 
terrain. As was previously pointed out, the Fortran code 
written for evaluation of line of sight in three dimensional 
terrain was not verified. With the availability of digital- 
ized terrain this portion of the program should be verified. 
However, one serious disadvantage in the concept used in 
this model is the requirement that the target course line 
be parallel to the X axis of the coordinate system. The 
problem lies in the necessity of assigning new altitude 
values for each different route of approach into the target. 
Another method for determining line of sight between 2 
points in 3 dimensional terrain is found in reference 6. 

This method requires a least square fit to the terrain 
points lying on the line between the two points for which 
line of sight is being determined. Then the method deter- 
mines if the curve extends above the line of sight. 
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V, 
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INATES OF MA 


X 


H 


0.0 


6066. 


6706. 


6123. 


9078. 


5757. 


19087. 


5426. 


30C00. 


5258. 


48854. 


5310. 


52714. 


5556. 


56446. 


5978. 


58837. 


6097. 


60886. 


6075. 


6220C. 


5890. 


65245. 


6048. 


68374. 


5819. 


72325. 


6374. 


78000. 


5*50. 


84023. 


5630. 


85762. 


5 504. 


92782. 


5793. 


97683. 


6096. 


98753. 


6342. 


100321. 


6494. 


101633. 


6389. 


10508C. 


6617. 


106 60 C. 
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7328. 
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8160. 
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7423. 
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251296. 


7068. 


257179. 


6792. 
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6919. 


262358. 


6750. 


264328. 


7021. 


265230. 


7460. 



JOR TERRAIN POINTS 
Y 

0.0 
0.0 
o.c 
0.0 
0.0 
0 . 0 
0.0 
0 . 0 
0 . 0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
c.c 
c.o 
0.0 
c.o 
0.0 
0.0 
0.0 
c.o 
o.c 
o.c 
0 . 0 
0.0 
o.c 

з. C 
0.0 
c.o 
0.0 
c.c 

и. 0 

0 . c 

O . 0 

c.o 
0.0 
0.0 
J.O 
0.0 
0.0 
0.0 
0.0 
0 . 0 
0.0 
•o.c 
0.0 
0 . 0 
o.c 
0.0 
c.o 
o.c 
o’. 0 
c.o 
o.c 

G • 0 
0.0 
0.0 
0.0 
c.o 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0 . 0 
, 0.0 
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COORDINATES OF MAJOR TERRAIN POINTS 
CONTINUED 



265927. 


7496. 


C.O 


268257. 


7028. 


0.0 


275105. 


6983. 


0.0 


279714. 


7203. 


0.0 


281700. 


7300. 


0.0 
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8045. 


0.0 


286963. 


8012. 


C.O 


m 61 : 


UlS: 


h:8 


291700. 


7162. 


c.c 


294290. 


7363. 


L.C 


295721. 


7773. 


C.O 


297609. 


7534. 


0.0 


300001. 


7600. 


0.0 


302107. 


7627. 


0.0 



COORDINATES OF 


X 


Y 


196703. 


0.0 


212183. 


0.0 


22440C. 


0.0 


265927. 


0.0 


159925. 


0.0 


18920C. 


0.0 


226653. 


0.0 


259793. 


0.0 


286963. 


0.0 


290001. 


0.0 


285209. 


0.0 



BATTERY POSITIONS 



795C. 

8160 . 

7496. 

7826. 

6700. 

8312. 

6919. 

8012. 

72t>5. 

8045. 
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APPENDIX III 



The 


GLOSSARY OP FORTRAN VARIABLE NAMES 
following is a list of Fortran variable names in 



the order of appearance. This list excludes the system 
parameters which were presented in Chapter III. 



K(K,LY) 


THE X COORDINATE OF THE KTH TERRAIN 
POINT WITH THE LY TH Y COORDINATE INDEX 


H(K,LY) 


THE H COORDINATE OR ALTITUDE OF THE KTH 
TERRAIN POINT WITH LY TH Y COORDINATE 
INDEX 


Y ( K , LY ) 


THE Y COORDINATE OF THE KTH TERRAIN POINT 
WITH LY TH Y COORDINATE INDEX, SET EQUAL 
TO 0 FOR TWO DIMENSIONAL TERRAIN 


BX(J) 


X COORDINATE OF THE JTH BATTERY 


BY( J ) 


Y COORDINATE OF THE JTH BATTERY SET EQUAL 
TO 0 FOR 2 DIMENSIONAL TERRAIN 


BH ( J ) 


H COORDINATE OF THE JTH BATTERY 


RANGE ( I ) 


AN ARRAY USED WITH SUBROUTINE UTPLOT 


AH ( I ) 


ALTITUDE ABOVE TERRAIN FOR THE ITH POINT 


DELTX 


THE DIFFERENCE IN X COORDINATES OF THE 
KTH TERRAIN POINT AND THE JTH BATTERY 


DELTY 


THE DIFFERENCE IN Y COORDINATES OF THE 
KTH TERRAIN POINT AND JTH BATTERY 


DELTH 


THE DIFFERENCE IN THE H COORDINATES OF 
THE KTH TERRAIN POINT AND THE JTH BATTERY 
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ELRX 

AZRX 

RMF 

RM 

THETA 

AT 

TR 

RS 

HF 

RXS 

HFT 

RSA 

RSS 



MAXIMUM ELEVATION TRACKING RATE FOR THE 
SPECIFIC BATTERY BEING EXAMINED 
MAXIMUM AZIMUTH TRACKING RATE OF THE 
SPECIFIC BATTERY BEING EXAMINED 
MINIMUM RANGE OF THE SPECIFIC BATTERY 
BEING EXAMINED 

MAXIMUM RANGE OF THE BATTERY BEING 
EXAMINED 

DETECTION RADAR LOWER EDGE OF BEAM 
PATTERN ELEVATION ANGLE IN DEGREES OF 
THE BATTERY BEING EXAMINED 
SYSTEM DELAY TIME FOR THE BATTERY BEING 
EXAMINED 

CREW REACTION TIME FOR THE BATTERY BEING 
EXAMINED 

SLANT RANGE TO THE (K,I)TH POINT 
HORIZONTAL COMPONENT OF MINIMUM RANGE 
HORIZONTAL COMPONENT OF RS 
HORIZONTAL COMPONENT OF THE RANGE OP THE 
POINT WHERE THE LINE DRAWN FROM THE 
MAXIMUM ELEVATION ANGLE INTERSECTS THE 
TARGET COURSE LINE 

THE SLANT RANGE WHERE MAXIMUM AZIMUTH 
RATE OCCURS 

THE SLANT RANGE WHERE MAXIMUM ELEVATION 
RATE OCCURS 
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XD 


THE DISTANCE ON THE TARGET COURSE LINE 
AVAILABLE FOR FIRING ALL BUT FIRST AND 
LAST ROUNDS 


HOP 


OPEN FORE SLANT RANGE 


T 


TIME AVAILABLE FOR FIRING EXCLUDING 
FIRST AND LAST ROUNDS 


NS ,XNR 


NUMBER OF ROUNDS FIRED PER RUN 


RH 


HORIZONTAL RANGE TO MAXIMUM RADAR 
DETECTION POINT ON COURSE LINE 


TP 


TIME OF FLIGHT 


GUNS, NGUNS 


NUMBER OF GUNS PER BATTERY, RADAR CON- 
TROLLED GUN BATTERY 


R(II) 


RANGE OF II TH ROUND FIRED 


SIG(II) 


THE FIRING ERROR RADIUS OF THE II TH 
ROUND FIRED 


SIGM 


FIRING ERROR OF THE RADAR CONTROLLED GUN 


EP(II) 

Ti 


THE PROBABILITY OF KILL OF THE II TH 
ROUND FIRED 
NS 

7T (l-EP(II)) FOR THE RADAR CONTROLLED 
L=I 

BATTERY 


PD(J) 


PROBABILITY OF KILL OF THE JTH BATTERY 


TY 


NS 

7T (l-EP(II)) FOR AA'rf BATTERY 
L=I 


NT, XNR 


NUMBER OF MISSILES FIRED 


TP 


NFU 

TT (l-PK(J)) THE PROBABILITY THAT NO 
D=I 

ROUND KILLS THE TARGET 
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PHI (K, I ) 


PROBABILITY THE AIRCRAFT IS KILLED AT THE 
( K , I ) TH POINT IN SPACE 
SUBROUTINE LOS 


RXM 


HORIZONTAL COMPONENT OF MAXIMUM FIRING 
RANGE 


XNK 


X COORDINATE ON TARGET COURSE LINE OF 
FIRST RADAR DETECTION POINT 


KK 


X COORDINATE INDEX OF LAST TERRAIN POINT 
BETWEEN POINT IN SPACE FOR CALCULATION 
AND THE BATTERY 


II 


X COORDINATE INDEX OF FIRST TERRAIN POINT 
BETWEEN POINT IN SPACE FOR CALCULATION 
AND THE BATTERY 


SLPD 


THE SLOPE OF THE LINE BETWEEN THE POINT 
OF FIRST RADAR DETECTION AND THE BATTERY 


SLP(M) 


THE SLOPE OF EACH OF THE M TERRAIN POINTS 
THAT MAY CAUSE TERRAIN MASKING ON EACH 
FIRING RUN 


Z 


DISTANCE ON TARGET COURSE LINE FROM POINT 
OF TERRAIN UNMASKING TO THE BATTERY 


XDR 


THE DISTANCE BETWEEN POINT OF FIRST 
RADAR DETECTION AND OPEN FIRE POINT, ALSO 
CALLED WARNING TIME DISTANCE 


TT 


TOTAL DELAY DISTANCE ON TARGET COURSE 
LINE SUBROUTINE SEE 
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EYE 



INDICATION OF WHETHER LINE OF SIGHT EXISTS 



FOR EACH CALLING OF THE SUBROUTINE SEE 
EYE = 0 INDICATES LINE OF SIGHT EXISTS 
EYE =1.0 INDICATES LINE OF SIGHT DOES 
NOT EXIST 
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APPENDIX IV 



FORTRAN PROGRAM 
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DATE * 63325 
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I f ( ( xrunp t T x ) , lt* { wi *vr* t i o » ) go r n n* 

p ( ! M =DSQPT ( (XOfOEl TX-WT*VT*T1P) «*o f nEl TH*DF( TH*nFI t y*OFI TV) 

SIG( I I ) = ( ,^^*VT)7 + 7f(p(Tl ) **2*1 n , OF -6 ) t ( $ I GM^VT/ 1 0-7 n ,A)7*n( 
i A P I + Y M* p ( I I ) ) ) * * ? 

FP(T!) = fA/(?.0^'1.nno?S c *STG{TT)))-PG 

TW=Tw*( 1.0-Frni ) ) 

P«=P ( I I ) 

S I G 1 = S ! G ( T ! ) 

C P! = FP( T T ) 

IRF rONTTNUF 

PK ( I ) = ( 1 . n-( TW» **NGUM) 
r ^IMT All, TW.J t PK( J) 

r 4ii wnrx.'li-nin-P?!, ,.u- n N) = ',r,i2 4,//, 
r n^x,* r* ( • f i 2 »•)=*, ix f pi".Q) 

0 DO TMT f X\'P ,T, T F 

r R01 r n9^U< MX, »NU*RFR OF ROUNDS P T® MX , F T . ^ , SX , • T = • , 1 X , r,K . ? t ^ , 

G l , TF- , ,lX,F1 r .4) 

GO T n ir 

r automatic w c £pon r^.ttpoy 

60° T =( XO/v T ) 

MC=TitTCR A 

IFfPH.GT.PYM) NS=NS*1 
f HMS.G T .S rt ) NS = «M 
XNP=NS 



r n s 

H H 

* m 
•m 



1 1 o 



'■I 70 
a 1 n 1 

' 1 n 

"i?R 
f 74 
~ 1 *> R 
26 



•1 27 
“ 1 2P 

A | po 

n 17 n 
"hj 
^1 ?? 
^ 1 *3 

^174 

M1F 

MIA 

'm 



r HR 
*<1 in 
MP 

: 147 
"1^3 
n 7 44 

a i 4 p 

' ! 4 6 

M47 

c ' 4R 
m 40 
MFn 

'in 

^ I c 2 

A 7 q 3 

- 1 c 4 

- 7 c*; 

MF7 
n 1 « P 



TY=1 . 0 
NAGUN=AGMN 

OP M44 (1=1, 

W I = T T - 7 

IF ( ( XO+OFL tx ) .LT. ( WT *VT* i . ^ /TCP \) ) Gn th 144 

D(TT)=ncnp^(( xn^OFt T X-WT ♦ VT* l . ^ / T CP A )** “> ♦ OFl th *• OFLTH^nFi TY*OEt TY) 
S TGM 1 ) =( . 'MVTtcQ (iinn, "«--R )-*•> 

F P ( f I |-(A/(?. r '*R, 141EP?65Y C 1G(1I M)*PA 
TY= TY* ( 1 , "-CP( IT)) 

I 44 rPNTIWJF 

° K ( n =( 1 . 1 -( TY n*N*GlJN» 
r OQ PJT 410 , ty f J f OK ( J) 

0 4n r F>pMAT( x , • ( 1 -on ( 1 -P^> ) . . . ( I -PM) *■ ,GIP- 4 , // , 
r mx f 'PV ( • * 12 f * ) = » f 1 K,FP, Q) 

r PPTNT 640 t X^ P f o<( )) 

C A4.~ PPR^AT{ 10 X , * XMO = * , F? •' , “ , 1 X , • 0K( .) ) =* , 1 x,f 70. 7 , PX , • A A PT F V • ) 

GO rn ■> a 

r MTSSTIF WOHFI 

T= ( X 0 / V T 1 
WP=MP 
NT = mp 

!F< T.GT. ( W0*T1M) ) GO To ip 
NS=T/T! m 

xnr=m<; 

TPMP = 1. n 

X^P = 0 M A XI { TC VP , XNP ) 

N T =XNP 

TP OK ( m 1 . P- ( 1 . O-FXSM) »*^JT 
r DO T NT A,''*', T t TF » c FO 

0 cpRM a t ( ’ r X , • T=- • f r 1 r . 7 . ♦ rr = • , F 7 - l f 'ccd = i ,r,p. ?) 

0 PPTNT t-c; f MP,r t PK| ]| 

f to c C0P M AT ( 7 r x , • NP = f , 1 X,n f lX* ,T =**'X f Gl '.2,Fx f tMTCcfLCt t o,'F t 4) 

GP jn R a 
1 0~ PK( n=n, 

7 * '"PNMNUe 

TPs], n-OY ( 7 ) 

Nl=N c UM 

7P 7 on JTa? t K‘L 

TP = t04c( 7 . n-pK( )] ) ) 

7 Q“ CP NT* I vnr 

PHMX,T ) = i,r-TP 
AXH( x t t ) =A m( t ) 

7^ CPNTINUF 
IT OONTTNitr 

OP 7RA f =nti , nt t 
OP 761 >= 1 ,^ r 

A A = 1 
ALT = AA< ' 

NLT = .U T 

PHI ( 11 = ( PH! ( ! , in(I. r -PO(J)|) 
ni GPNT INIJF 
P 0 I N T T S 2 
7S7 FPP M A T ( 1 H 1 ) 
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0 16Q 
n 16" 
0161 
"16? 
0163 
01*4 

0 166 

01 66 
"1 67 
"168 
" 1 69 
r 1 *T 7 

6 1 “M 
cnpTPA^j 

onr.i 



r r *3 



'f'C4 
> r» r « 
apt * 
•'• *• r "» 
"TR 

AP ~ Q 



* " i r\ 

'*>11 
P 1 7 

vi ^ 

-A , ^ 
-' r 1 6 

' r ’ 6 
*\ r ^ -r 

.rip 
** " * o 
?r 
'.a 7 i 



"":>3 

A A A 4, 



f "?6 
AP ?* 
A A p 7 
r a 7 p 

a a p Q 



•" a 



v “ A 7 » 

'"■n 

a 7 7 
AA 74 

A A O ^ 

r a 7* 

A *% -j 7 
.' A 7 p 

Ar 30 
A A4.A 
“A41 

r 047 

r A 4 3 
aa 44 

; a 4 6 
P O 4 4 

A A 4 7 
r a 4 p 
A a 40 

PA 

'‘AC) 
A A c 7 
A A*p 
r a c 4 

A a t; q 

AA^^ 

0 A c 7 
•* P * 0 
A 7P p 
0 a *n 

p a ^ ■. 
AA* 7 
A A *7 

m 64 

PA 4 q 

r "66 
r a A 7 

706 3 

A A 4 O 



Ap 7P 
A P 7 t 
A p 7 7 
CP77 



75" CONTINUF 

00 176 K=NT1 f NT- 
on 175 ! *1 ,30 
AA»T 



4X1 ( I ) = AA*10O. 

PHI ( 1 ) = PHIIK f n 

176 CONTINPP 
PRINT 177. K 

177 F0PMAK 1H1.67X, MLTUTinF VS PPO R » , / / , 64 X 
CALL 1JTPL0T( AXI,PH1 ,7",RANGe t i f n) 

176 CONTINUE ? 

Ifil CONTINUE 
FNO 



’TPRPAIN PT« , 7X, 13, //) 



lv n t F\'Fi 1 1 mod 7 ins 



n AT F = 63 3 TP 1 7/61 /66 



r 

f 



r 

r 

r. 

r 

r 



r 

r 



r 

r 



r 

r 

r 

r 

C 

r 



r 

r 



6 

1 67 



3? 



71 

76 

7 

1° 

? A 

M" 



21 

p 

?0o 

1 43 
1 5 C 



74A 

35^ 

356 

727 

"2" 



5 o r 
60^ 
70- 

77 



56 

236 

237 

20 

27 

323 



SURPHUMNF LnS(PM,*.l t ) , l Y ,0F t TX , OF ). T Y , HP LTM , AO I , YM f T HF T A , TP , AT , V T 
* ? YK f XO. TF f TX, NTT .LT1 .LT? f I. T* ,f "NT f 7M f SM f PH.PXM, RS* > 

Implicit R PA ) * 3 ( A -H ,0-7 ) , INTFGrp*4( T -N) 

COMMON UH 3A») f H(1 70,u , X( I ?A f 1 ) f V( 17P , 1) f pxr ?n 1 ,3V( 7C 1 ,RH( 77 ) , p 

1 K ( 30 ) ,PMJ ( J.?0 30) .AXHd^n,^ ) ,6) P (OP) |0( 2 pr , ,SK,( ?aa ) , pp( ? pp , 

OAlCMLATF OT^TANCE ON COHftSF LINE cn » FIojmo 
C At C ML A T c RADAR WARNING T T M c 
BKs, ?1^ 

7K-. 24 

RXM = nS0RT ( PMtRM-nF) TH*DFL TH-DFt T V *n FL T Y ) 

0 H=OS0R T ( ( OFITH/DTANI THFT4) ) 7*7-0 FL T Y** ? ) 

XNK = RX(.n-n6GRT( RH*RH-O p LTV*nE( TV ) 

IF ( XNK,, LT, n. | XNK*C . 

PRINT 6 , K , I f J 

F0R*4T( 10X,'TP=',1X,T3,TX,'ALT=«,2X,I3,1X,'RTRV=',1X,I3) 
op I NT 153, XNK ,RX( J) ,PXM 

F0OM A T ( 10X, 'XNK- • ,lX,Fl r '. n *6X, , RX(J) = , ,lX,F1c ,'l t l X, • PX M=I ?1X,F1"#C 
1 ) 

f c( xnk. fo. ° x r j ) ) on to 23 

NX K = XNK 

T F ( R Y ( J I.NC.a. ) co T n 
11=0 

no 3i Ms ] ,ntt 

I F ( I I.NF. p) c.n TCI 7 2 

I F ( XNK, LT, X( M .L V ) ) I I = M 

I F ( B X ( J). C,F. X? M f LY11 GO TO 71 

KK=M-7 

GT Tn 3 6 
CONTINIIF 

SL PF = ( AH ( I 1 - PH ( J) ) / ( R X ( J ) - X NK ) 

^P I NT 7 q.PP , K K , I I 

format ( iox,'SLP6=* t ix,Fir.5,lX, * KK = ',I3,MT = ',m 

DO 700 M=I I f KK 

I F ( KK- ! I ) 77 f i o ,10 

St P(M) = ( AH( n-HC M f ( Y n / ( X( M,LY) -XNK ) 

PRINT m, $LP( M ) » X ( M , L Y ) 

FOP M AT(!CX, , Sl P( ’ •I3, , l = , ,Fl r .R,lX,*X( M ,LY) = * , p 1 7 , ) 
IF(SLPR.GT.SLP( M ) ) GO TO 310 
GO TO 7pr 

Z = (AH(!)-RH(J))MRX( )) -X ( ^ , L Y ) ) / ( H t « , LV ) - RH ( I ) ) 

XNK =R X ( J) -? 

SI PB=( AH( I )-PH(J))/t^X( M -XNK ) 

PRINT ?l . XNK ,RX( )) , M 

FOR w A T ( icx, * XNK = * f F 1 r # 0, * RX( J) = * ,F] :7. n f *m = . , I 3) 

PRINT R,SLPP*AH(T) 

FOOMATt IPX, * S L PR=« ,F! C. 6 , ' AH( !) = • ,F10.r , 

CONTTNMF 
PRINT 143, XNK 

FORMAT ( !C X, • X NK = * , IX, GIG. 4 1 
GO TO 

KK = ( X(K,LY)-XNK) /YK-H , 7 
xmk=p. 

MNxO. 

n=". 

nn 320 M = ? , k k 

CALL $FF(nxk,k # i,j,eyF,LY,yk) 

I F f c Y F ) 32?,?2 7 ,74 r 

1F( MN.FQ. ( M -1 ))GH TT 360 
GO TO 766 
T S=71/VT 

lF(TS«G’ r . TX) XNK = XNK+ YK*XM< 

71=71 ♦ YK 
MN = M 



NXK=XNk 



XMK=xmk+ 1 . 0 
CONT T NM C 

oH=nSOPTf ( RX( J» -XNK ) *♦2-051 TV*DFLTY ) 

R X M = p m i N1 ( RXM t pm) 

Q0C=r)60R T ( RX^^RXM-fDCLTH^npLTH^nFLTY *OFL T Y ) 
7(1) =pnc 

IF(J.LF.LTl) GO tp 6AP 
I c ( J.LF.t T 2 ) GO TO 6^0 
1F( J.LF.I T^» GO Tp 7P0 

TPsAO (4-yMtPPP 

r,n to pc 

TF=OEXP(C0NTTRnc , 

GO TO 7 p 

TF= ( P OF- 7 * ) /SH 

T T= ( TR4.4i4.Tp ) ♦v T 

XOR=RH-PXM 



I f ( xnR, c 0 . 0 , ) r,n to 23 s 

1 F ( TT ,L T, xnp ) GP TP 66 
TT= TT-XOD 
GO TO 7 3 p 

xn = p XM-OFL TX-TF-k VT 
GO TO 737 

XD=R XM-DFLTX-TT-tf^VT 
C ON T T NM C 

PRINT ?n 9 OH,RXM,XD,XnR 
FOR M A T ( IPX, 

GO TO 373 

xn = 0 . 



,tt 

6G10. 4) 



RFTMRN 

r N o 



104 



FORTRAN IV G L p VEL MOD 2 



SEF 



DATE = 68X28 



0^01 
no^ * 

r\Ao > 

' X04 
opos 

^07 
f ' onq 

A a q 

"»o i n 
M 
A n\2 

''-X 1 X 

>r i 4 
5 

x* ^ * 

'AH 

AA ! 8 
TOI Q 



4" 



^UROnuTTNJF SFE<NXK t K.I .J f F yF.I Y.YK) 

Implicit REAl*8< A-H f O-£ > , TNTEGFR*4H-N> 

r DM yon AH( 3">) .H( 7?0,n ,xn 2D 1 1 ) , Y< 1?7 , 1) ,PX( ?x ) , RY(2 A I 
IK ( x- ) t PHI | 1?C , XO } f AXH< 12° t ^0 ) tSLPlO'M ,R ( ?0 r I ,8 IG( | f FP< 
XNK=NXK 

RMRx»DS 0R T ( ( PY( .1) -Yf K . LY M **2* I 8X ( J ) -XNK I *♦ ? > 

SL OR * ( H ( J 1 1 Y ) -RH ( J ) )/RMRX 
N*C(PY(J)-Y(K # LY) )/YK)-. 9 999 
NYK=YK 
L YK=L Y 

SS = ( RY< J ) -L YK ) /R MR X 
ORTMT 40. 

FCPMATMAX, * PTRYlS* f I3 f lox f FIO 4) 

DO n I I = 1 . N 

XK=YK*pT6N(nA RCO^ ( ( PY( J) -LYK ) /PMPX M 
L XKs ( XK + XNKI /YK 
y y wr — i yy 
OX*XK + XNJK-XXY 
IFCDX.l T . .*) GO TO 5 a 
L X k = L X* ♦ 1 



r\r\pr\ 


50 


L YK=L YK+ YK 


ro27 




X T - I ♦ 1 


A A 2 2 




SLPJ-(H( K f l Y ) - H ( LXK f LYK) ) / DS OR T (YKtYK*XK*XK ) 


~«xx 




T F f ci PP„ GT. M PJ ) GO TO ' r 0 


. 3?4 




YK - YK *X T 




l A 


COMT I nue 


~ * 2 fi> 




FYC=0. 


■‘027 




GO rn x>> 


^ \ ■> ft 


1 


c Y r = 1 , n 




V 


GO JTTMlic 


■\n^n 




«F T ORN 


-SA q | 




F NO 



1 7/SI /SS 
.RH( ,P 

y n I 
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APPENDIX V 



PROBABILITY VERSUS ALTITUDE EOR SELECTED 
TERRAIN POINTS 
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ALTITUDE VS PRDB 
TERRAIN PT U 




CURVE 1 



ALTITUDE VS PROB 
• TERRAIN FT 1»7 




CURVE 2 
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ALTITUDE VS PfiOB 
' TERRAIN PT U8 ' 




CURVE 3 
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ALTITUDE VS PROB 
TERRAIN PT 62 ' 




CURVE 4 



ALTITUDE VS BROB 
TERRAIN PT 69 




CURVE 5 



111 



ALTITUDE VS IROB 
. TERRAIN PT 82 




CURVE 6 



112 



ALTITUDE VS PRC© 
• TERRAIN PT 72 




CURVE 7 



ALTITUDE VS PBOB 
• TERRAIN PT 68 




CURVE 8 



APPENDIX VI 



COST VERSUS ALTITUDE FOR SELECTED 
TERRAIN POINTS 



115 



ALTITUDE VS COST 
TERRAIN PT Ul 



2.0 



1.6 



1.3 



1.0 



.67 



.33 




0.0 



750 



1500 



2250 



3000 



CURVE 1 
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ALTITUDE VS COST 
TERRAIN PT hi 




CURVE 2 



AUirma vs cost 

SBUIX PT W 




CUBVE 3 
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ALTITUDE VS COST 
TERRAIN PT 62 




0.0 750 1500 2250 3®0 



CURVE k 



U9 



ALTITUDE VS COST 
TERRAIN PT 69 




0.0 750 1500 2250 3000 



CURVE 5 



120 



ALTITUDE VS COST 
TERRAIN PT 82 




CURVE 6 
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APPENDIX VII 



DATA POINTS FOR THE CURVES IN APPENDICES 

V AND VI 
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TERRAIN POINT 41 



PROS 

0.0027158 
0. 002711C 
0.0954828 
0.4619637 
0 . 474601 J 
0. 5874727 
0.7439264 
0. 7306477 
0. 71 J6845 
0. 824807C 
0.9328496 
0.9329319 
0. 9559594 
C. 9597607 
0.9563827 
0.9864236 
0.9782999 
0.9858453 
0.9815221 
0.9825184 
0. 9809653 
0. 979332C 
0.9776287 
0.9758614 
0. 974C362 
0.9721588 
0.970235C 
0. 9683CC2 
0.9663471 
0.9643678 



ALTITUDE 
100 
2 00 
300 
400 
500 
60C 
700 
800 
900 
10C0 
1100 
120C 
1300 
1400 
1500 
160C 
1700 
1800 
1 90C 
2000 
2100 
2200 
2300 
2400 
2500 
2600 
2700 
2800 
2900 
3C0C 



COST 

0.6375157 
0.4982110 
0.4660828 
0.7422637 
0.7152010 
0.7920726 
0.9031264 
0.8490477 
0.8013845 
0.8931070 
C. 9937496 

0. 9874319 
1.0060587 
1.0072603 
1.0017824 
1.0305233 

1. CC90990 
1.0098448 
1.0030212 
1.0C35181 
0.9981653 
0.9940320 
0.9920237 
0.9901614 
0.9847361 
0.9792588 
0.9763350 
0.9736002 
0.9704471 
0.9676678 
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TERRAIN POINT 47 



PROB 


ALTITUDE 


COST 


0.0 


ICC 


0.6348000 


0.0 


20C 


0.4955000 


0.0 


300 


0.3706000 


0.2796825 


400 


0.5599825 


0.004294C 


500 


0.2448940 


0.0027300 


600 


0.2073299 


C. 005534* 


700 


0.1647344 


0.0027246 


80C 


0.1211246 


0.0027211 


9C0 


0.0934210 


0. 00271 7C 


1000 


0.0710170 


0.0027124 


1100 


0.0636123 


0.0027072 


1200 


G. 0572072 


0.0054274 


1300 


0.0555274 


0. 54 50826 


1400 


0.5925826 


0.5536741 


1500 


0.5990741 


0.393517C 


1600 


0. 4376170 


0.8372307 


170C 


0.8680307 


0.7275365 


1800 


0.7515365 


0. 8045451 


19C0 


C.826C451 


0.8778201 


2000 


0.8988200 


0.9350527 


2 ICC. 


0. 9522527 


0.9264631 


2200 


0.9411631 


0 .9693650 


2300 


0.9837649 


0.9780607 


2400 


0.9923607 


0.9817001 


2500 


0.9924001 


0.9815891 


2600 


0.9886891 


C. 9841095 


2700 


0.9902095 


0. 979245C 


2 800 


0. 984545C 


0. 9829284 


2900 


0.9870284 


0.9814092 


3000 


0.9847092 



124 



TERRAIN POINT 48 



PROB 

0.0027316 
0. C027297 
0.0045465 
0.0027241 
0.0027205 
0 • 00 2 71 64 
0.0027117 
0.4124928 
0.4975135 
0.5125919 
0.6919966 
0.7918245 
0. 801255C 
0. 784154C 
0. 8661683 
0.9242823 
0. 9468088 
C. 9326292 
0. 9726201 
0.9639952 
0.9749367 
0. 98*+32G6 
0.9811772 
0.98*2251 
0.9827261 
0.981188C 
0.9795651 
C. 9778711 
0.9761124 
0.974295C 



ALTITUDE 

100 

zee 

30 C 
400 
500 
600 
700 
800 
900 
1000 
1100 
1200 
1300 
1400 
1500 
1600 
1700 
1800 
1900 
2CCC 
2100 
2200 
2300 
2400 
2500 
26C0 
2700 
2800 
29CC 
3CC0 



COST 

0.6375316 
0.4982296 
0.3751464 
0.2830241 
0.24332C5 
0.2073163 
0.1619117 
0.5308927 
0.5882134 
0.5808919 
0.7523966 
0.8463245 
0.8513550 
C. 8316540 
0.9115682 
0.9683923 
0.9776088 
0.9566292 
C. 9941202 
0.9899952 
0.9921367 
0.9990206 
C. 9955771 
0.9985251 
0.9934261 
0.9882830 
0.9856650 
0.9831711 
0.9802124 
0.9775949 
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TERRAIN POINT 62 



PROB 


ALTITUDE 


COST 


0.0027238 


IOC 


0.6375238 


0.0027202 


200 


0.4982201 


0.0027160 


300 


0.3733160 


0.0027113 


400 


0.2830112 


0.4737283 


500 


0.7143282 


0.4753606 


600 


0.6799605 


0.4886442 


700 


0.6478442 


0.6693962 


800 


0.7877962 


0. 7629707 


900 


0. 8536707 


0. 7530464 


1000 


0.8213464 


0.7393761 


1100 


0.80C2760 


0. 8401864 


1200 


0.8946863 


0.90718G8 


1300 


0.9572808 


0.9384188 


1400 


0.9859188 


0.9611555 


1500 


1 . OC 65546 


0.9649823 


1600 


1.0090818 


0.9611157 


1700 


0.9919156 


C. 9878440 


180C 


1.0118437 


0, 980627C 


1900 


1.0021267 


0. 9793023 


2CC0 


1.0003023 


0.9825637 


2100 


0.9997637 


0.9825942 


2200 


0.9972941 


0.981C45C 


2300 


C. 9954450 


0.9794154 


2400 


0.9937154 


0.9777154 


2500 


0.9884154 


0.9759512 


2600 


0.983051! 


0.9741287 


2700 


0.9802237 


0.9722538 


2800 


0.9775538 


0.9703322 


290C 


0.9744322 


0. 9663934 


3000 


0.9716933 
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TERRAIN POINT 68 



PROP 

0.1368645 
0.002731 5 
0.C027295 
0.0G41496 
0.0C27739 
C.0C27203 
0.0027161 
G. 0027114 
0.4335030 
0.4805196 
C.49413C2 
0.6767096 
0. 77CG04A 
0.762040? 
0. 7501766 
0. 8461959 
0.91’ 2051 
0.9404540 
0. 9631043 
0.9668788 
0.9629577 
0.9883308 
C. 9815130 
0. 9797094 
0.9829618 
0. 982624* 
0. 9810769 
0.9794487 
0. 9777501 
0. 975987C 



ALTITUDE 
ICO 
200 
300 
400 
500 
600 
700 
80C 
90C 
1000 
1100 
1200 
1300 
1400 
1 50C 
1600 
1700 
1800 
1 90C 
2000 
2100 
2200 
2300 
2400 
2500 
2600 
2700 
280C 
2900 
300C 



0.7716645 
C .4982314 
0.3733295 
0.2844496 
C. 2433239 
0.2073202 
0.1619161 
0.1211113 
C. 5242029 
0.5488196 
C. 555C302 
0.7312096 
0. 82C1044 
0.8095402 
0.7955766 
0.8902959 
0.9420050 
0.9644540 
C • 9846043 
0.9878788 
C. 9801577 
1.0030308 
0.9959130 
0.9940094 
0.9936618 
0.9897244 
0.9871768 
0.9847487 
0.9818500 
C. 9792870 



12 ? 



TERRAIN POINT 69 



PRCB 

0*0027225 
0. 002 71 87 
0.0027143 
0.0027094 
0.0066825 
0.3684522 
0.3379564 
0.8C41257 
0.7355692 
0. 8131848 
0.7832339 
0.8648104 
0.8990118 
0.9483667 
0.9563300 
0. 9776988 
0.9644131 
C. 9816354 
0.9819365 
0.9795866 
0. 9827996 
0.9820328 
0.9804507 
0.9787943 
0. 9770698 
0.9752833 
0.9734408 
0. 97i 5480 
0.9696105 
0.9676816 



ALTITUDE 
100 
200 
300 
400 
500 
600 
70C 
800 
900 
1000 
1100 
1200 
1300 
1400 
1500 
1600 
1700 
1 80C 
1900 
2000 
2100 
2200 
2300 
2400 
250C 
2600 
2700 
2800 
2900 
3000 



COST 

0 . 
C. 
0 . 
0 . 
0 . 
0 . 
0 . 
0 . 
0 . 
0 . 
0 . 
0 . 
c. 
0 . 
1 . 
1 . 
0 . 
1 . 
1 . 
1 . 
0 . 
0 . 
0 . 
0 . 
0 . 
0 . 
0 . 
0 . 
0 . 
0 . 



6375225 

4982187 

3733143 

2830094 

2472824 

5730522 

4971564 

9225257 

8262691 

8814847 

8441339 

9193104 

9491118 

9958667 

0017300 

0217981 

995213C 

0056353 

0034361 

0005856 

9999996 

9967327 

9948506 

9930943 

9877698 

9823833 

97956-08 

9768479 

9737105 

97C9816 
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TERRAIN POINT A? 



PROB 


ALTITUDE 


COST 


0.0027197 


100 


0.6375197 


0.C027155 


200 


0.4982154 


0. 00271C7 


300 


0.3733107 


0.0316351 


40C 


0.3119351 


0.445118C 


500 


0.6857180 


0.4572794 


600 


0.6618794 


0. 5182828 


700 


0. 6774828 


0. 7186771 


800 


0.837C770 


0.7030515 


900 


0.7937515 


0.6739765 


100 c 


0.7422764 


0. 8060656 


1100 


C. 8669655 


0.9227645 


1200 


0.9772644 


0.9237672 


1300 


C. 9738672 


0.9685816 


1400 


1.0160909 


0. 9673398 


1500 


1.0127392 


C. 9683690 


1600 


1.0124633 


0.9843186 


1700 


1.0151186 


0.9753590 


1800 


0.999359C 


0. 9849095 


1900 


1.C064C37 


0.9802112 


2000 


1.0012112 


0.9824312 


2100 


C. 9996311 


0. 9808694 


220C 


0.9955694 


0.9792317 


2300 


0.9936317 


0. 9775244 


2400 


0.9918243 


0.9757534 


2500 


0.9864534 


0.9739249 


2600 


0.9810249 


0.9720447 


2700 


0.9781446 


0.9701182 


2800 


0.9754182 


0.9681811 


2900 


0.9722811 


0.9662319 


3000 


0.9695319 
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